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Introduction 


The  structure  and  activity  of  the  solar  atmosphere  are 
dominated  by  magnetic  fields.  These  fields,  arising  from  sources 
beneath  the  photosphere,  channel  the  solar  wind  in  coronal  holes 
and  confine  the  high-pressure  plasma  in  active  regions. 

Of  particular  interest  is  the  ability  of  stressed  magnetic 
fields,  usually  with  both  ends  (foot  points)  embedded  in  opposite 
polarity  regions  on  the  solar  surface,  to  store  significant 
amounts  of  energy.  The  ground  state  of  a  magnetic  configuration 
coexisting  with  a  highly  conducting  plasma  is  one  in  which  no 
plasma  currents  are  flowing.  When  currents  are  driven  through 
the  field/plasma  structure,  due  to  relative  motions  of  the  field 
sources  or  to  voltage  differences  between  the  foot  points,  the 
magnetic  energy  increases. 

The  stored  magnetic  energy  in  active  regions  comprises  the 
preflare  state  (Van  Hoven  et  al .  1980) .  The  excess  (non¬ 
potential)  energy  in  this  state  must  be  sufficient  to  provide  the 
energetic-particle,  flow,  and  thermal  outputs  of  the  flare. 

Using  numerical  codes,  we  compute  the  structure  and  excess  energy 
content  of  the  coronal  magnetic  field  in  the  preflare  state  when 
the  current  is  generated  by  a  twisting  of  the  field  foot-points. 

As  stated  in  the  proposal,  the  objectives  we  hoped  to 
accomplish  with  this  funding  were:  (a)  "to  complete  development 
and  check-out  of  the  numerical  code"  and  (b)  "to  conduct,  at 
least,  preliminary  studies  of  the  magnetic  energy  build-up  due  to 
photospheric  motion."  We  were  successful  in  these  efforts  and, 
in  the  remainder  of  this  report,  will  discuss  some  specifics  of 
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code  development  and  the  preliminary  results. 

Model 

We  simulate  the  field  evolution  of  the  magnetic  loop  using 
cylindrical  (r,0,z)  geometry.  This  is  a  reasonable  approximation 
for  a  narrow  loop  (Van  Hoven  1981).  We  take  all  variables  to  be 
independent  of  the  0-direction  so  that  the  solution  is  confined 
to  the  r-z  plane.  We  further  assume  that  the  atmosphere  is 
initially  in  static  equilibrium  with  an  embedded,  axial,  current- 
free,  magnetic  field,  that  gravitational  forces  can  be  neglected 
when  compared  to  the  dominant  magnetic  forces,  and  that  the 
density  does  not  change  with  time.  Gravity  would  obviously 
influence  the  evolution  of  the  field  through  its  effect  on  the 
dynamics.  Compressibility  should  be  relatively  unimportant,  at 
least  during  the  early  stages  when  the  evolution  is  comparatively 
slow.  Both  of  these  quantities  can  be  included  at  a  later  phase 
of  this  study.  We  do  not  believe,  however,  that  either  gravity 
or  compressibility  will  significantly  modify  the  basic  phenomena 
we  are  investigating;  namely,  the  build-up  of  magnetic  energy  and 
shear  and  the  possibility  of  field  reversal.  The  plasma  0  is 
assumed  to  be  small  enough  that  the  pressure-gradient  forces  can 
be  neglected.  For  this  model,  the  single-fluid  MHD  equations  to 
be  solved  can  be  written  in  cgs  units  as 


along  with  the  auxiliary  equations 
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X  S, 


V  ‘V 
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The  total  derivative  is  represented  by  d/dt,  v  is  the  fluid 
velocity,  3  the  magnetic  field,  and  3  the  current  density.  The 
constants  are  the  resistivity  r),  the  speed  of  light  c,  the 
density  p,  and  the  viscous  coefficient  v.  The  viscosity  is 
introduced  for  numerical  stability  reasons  and  to  accelerate  the 
relaxation  to  a  static,  non-potential  equilibrium  configuration. 

The  magnetic  field  is  written  as 


5  .  i 


where  is  the  initial,  spatially  constant,  axial  field.  The 
radial  and  axial  components  are  then  given  in  terms  of  the 
magnetic  potential  'i'(r,z);  i.e.. 
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1  ^ 
r  Sz ' 


1  ^ 
r  dr 


Equations  (1)  and  (2)  are  made  dimensionless  in  terms  of  a 
distance  _a^,  the  hydromagnet ic  time  and  velocity  the 

initial  magnetic  field  and  the  density  p.  The  distance  a  is 
taken  to  be  the  total  radial  extent  of  the  photospheric  motion  in 
the  model  (this  will  become  clear  later).  The  three  time  scales 
in  the  problem  are  the  hydromagnet ic ,  resistive,  and  viscous 
t  i.res  g  i  ven  by 


These  time  scales  give  rise  to  two  important  parameters;  the 


Lundquist  number 
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and  a  viscous  number 
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where  a  simplified  form  of 


i  (r  * 

r  dr  dr^ 


dz‘ 


has  been  used  in  (3)  -  (6)  to  provide  dissipative  relaxation. 
The  advection  terms  have  been  dropped  in  view  of  the  small 
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velocities  (~  10”  a/ t^)  applied. 

Solution  Procedure 

The  method  of  excitation  used  to  simulate  the  photospheric 
input  is  to  rotate  the  ends  of  a  cylindrical  loop  (contained 
within  the  initial  axial  field)  at  a  rad ial ly-dependent  rate  out 
to  som.p  maximum  radius  ^  beyond  which  the  ends  remain  fixed.  The 
time  dependence  of  the  rotation  can  have  any  specified  form; 
i.e.,  the  rotation  can  be  increased  with  time  to  some  maximum  and 
then  either  remain  constant  or  be  decreased  to  zero;  the  rotation 
can  remain  fixed  at  a  given  value  throughout  the  computation;  or, 
alternatively,  it  can  be  reduced  to  zero  after  a  given  time 
interval,  etc.  Throughout  these  preliminary  computations  the 
rotation  rate  has  been  held  fixed  at  a  specified  level  during  the 
simulation. 

The  radial  variation  of  the  input  rotation  that  we  have  used 
to  obtain  the  results  presented  herein  is  as  shown  on  the  left 
edge  of  Figure  1.  The  azimuthal  velocity  increases  linearly  with 
radius  (as  in  rigid  rotation)  from  zero  at  the  cylinder  center  to 
a  maximum  (10  km  s"^)  at  r  =  0.8a,  after  which  it  linearly 
decreases  to  zero  at  r  =  a. 

The  remaining  quantities  to  be  specified  are  the  boundary 
conditions.  These  are  trivial  at  the  loop  axis.  Linear 
extrapolation  is  used  at  the  maximum  radius,  although  this  is  so 
large  the  variables  do  not  undergo  much  change  there.  For 
simplicity  in  these  initial  computations,  the  ends  of  the  loop 
are  assumed  to  be  twisted  identically  in  opposite  directions  so 


the  axial  mid-point  becomes  a  plane  of  symmetry.  Consequently, 
at  the  the  mid-point  (considered  to  simulate  the  loop  top), 
'?,Bq,u,  and  V  are  symmetric  and  w  =  0.  At  the  loop  base,  v  and 
Bg  are  given  by  the  specified  rotation,  u  =  ¥  =  0,  and  w  is 
linearly  extrapolated. 

Equations  (3)  -  (7)  are  now  solved  numerically  by  first 
solving  equation  (6)  implicitly  for  ¥  using  velocities  from  the 
previous  time  step.  The  remaining  equations  are  then  solved 
using  the  leapfrog  Dufor t-Frankel  scheme  (Peyret  and  Taylor 
1983).  The  leapfrog  method  was  initially  used  on  all  equations, 
but  numerical  instabilities  developed  in  ¥. 

Numerical  Results 

The  physical  values  used  for  the  results  presented  below  are 
as  follows:  a  =  10^  cm,  B  =  10.5  G,  p  =  n  m„  with  m,.,  =  proton 

o  '^CDOOp 

mass  and  n  =  10^  cm“^,  and  T  =  1.6x10^  °K.  These  parameters 
provide  p  =  0.1.  The  Lundquist  and  viscous  numbers  are  taken  to 
be  10^.  Some  features  of  the  solution  after  12,000  time  cycles 
or  approximately  260  hydromagnet ic  times  are  presented  in  Figures 
1,  2,  and  3. 

Three  field  lines  are  shown  in  Figure  1  as  they  would  appear 
in  projection  onto  the  r-z  plane  when  the  cylinder  is  viewed  from 
the  side.  Two  of  the  field  lines  start  at  the  base  within  the 
region  of  rotation  (r  <  a)  and  one  (the  short  dashed  curve) 
starts  beyond  it.  All  three  are  initiated  at  the  base  in  the 
plot  on  a  vertical  axis  (from  the  illustrated  viewpoint)  taken  to 
be  at  0  =  0°.  The  total  rotation  or  twist  of  each  of  the  field 


lines  be*-wee!i  the  base  and  mid-point  is  indicated  on  tlu^ 

The  field  line  starting  near  the  maximum  rotation  rate  at 
the  base  naturally  has  the  largest  amount  of  twist. 

T  n  r  e  r  s  t  i  ng  1  y ,  however,  more  than  half  of  the  twist  for  t  iu,'  field 
line  (122'^)  occurs  within  appr  ox  i  ma  t  e  1  y  0.2  d  i  imni;;  i  on  1  es.s  nr)  i  »  s 
(z  ',))  of  the  base.  The  inner  field  line  has  very  1  i  t  t  1  twi'd  , 
wfi  i  1  e  the?  outer  one  rotates  almost  uniformly  .iloivi  its  enl  it.' 

1  eng  th . 

The  radial  variation  of  the  magnetic  field  components  are 
plotted  in  Figure  2  at  the  axial  location  (z  =  0.8a)  indicated  on 
Figure  1.  The  most  interesting  feature  here  is  the  concentration 
of  the  axial  field  near  the  loop  axis.  This  is  seen  in  all 
simulations,  for  various  parameters,  that  have  been  conduf:t<'d  to 
date.  The  radial  component  has  a  maximum  near  the  axis,  .md  tin.' 
aziumuthal  component  peaks  somewhat  further  out. 

The  axial  run  of  the  field  components  are  shown  in  Figure  3 
at  a  radial  location  (r  »  0.8a)  near  the  region  of  maximum 
rotation  at  the  base  (see  Figure  1).  All  components  are  largest 
near  the  base.  The  axial  currents  generated  by  the  base  rotation 
feud  to  keep  the  .i/.imuthal  field  almost  const, mt  along  the 
cylinder. 

A  few  field  lines  from  an  analytic  solution  (Zweibel  and 
Boozer  1985)  to  a  problem  which  is  similar  to  fh,it  we  tiave 
considered  are  given  in  Figure  4.  Only  the  radial  and  axial 
field  components  are  shown  in  this  presentation  in  contrast  to 
the  remilfs  in  Fioure  1,  in  wh  i  c-h  t  ht*  aziinuthal  component  was 
also  consi.)lered  in  order  to  give  the  true  projection  view. 


7 


Although  Zweibel  and  Boozer  (1985)  had  to  make  some  additional 
approximations  (e.g.,  long  and  narrow  cylinder)  to  obtain  the 
analytic  result,  it  is  clear  that  our  results  and  theirs  are 
qualitatively  similar.  Both  show  that  most  of  the  radial 
curvature  occurs  near  the  base  (left  edge  in  Figure  4),  and  both 
exhibit  a  concentration  of  field  near  the  cylinder  axis. 

The  initial  magnetic  energy  contained  within  the  cylinder  in 
which  the  solution  was  computed  was  1. 11x10^^  ergs.  The  energy 
when  the  solution  was  terminated  was  1.29x10^^^  ergs,  representing 
and  increase  of  more  than  an  order  of  magnitude. 

Summa  r  y 

This  rese.irch  iias  demonstrated  that  tViere  arc  three 
character i St ic  features  which  occur  in  initially  axial  magnetic 
fields  as  a  result  of  twist  at  the  end  plates  of  a  cylindrical 
portion  of  this  field.  These  appear  to  be  very  general 
characteristics  since  they  occurred  to  a  varying  extent  in  all 
simulations  conducted  to  date.  First  of  all,  most  of  the  twist 
occurs  near  the  cylinder  base.  The  rotating  end  plates  also  tend 
to  concentrate  the  magnetic  field  near  the  cylinder  axis  over 
most  of  the  tube.  Finally,  the  largest  radial  curvature 
concentrates  near  the  base  in  conjunction  with  the  inward 
jiinctiing  of  tiie  field  lines. 

Fermi  ts  of  this  ctjsearch  program  have  been  presented  at  the 
annual  meeting  of  the  Division  of  Plasma  Physics  of  the  American 
Pliy;;  1  c.i)  Society  (Steinolfson  and  Schnack  1985).  We  also  intend 
t';i  present  some  of  this  work  at  the  summer  meeting  of  the 
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Figure  Captions 


Figure  1. 


Figure  2. 


Figure  3. 


Figure  4. 


The  rotation  velocity  is  sketched  on  the  left.  Three 
field  lines  are  shown  as  they  would  appear  in  projec¬ 
tion  when  the  cylinder  is  viewed  from  the  side.  All 
of  the  field  lines  start  in  a  vertical  plane 

(0  =  0°)  at  the  base.  The  inner  field  line  (— - ) 

has  very  little  twist.  The  center  line  ( - )  has 

the  most  twist,  but  the  majority  of  it  occurs  near  the 
base.  The  outer  field  line  {----)  is  twisted 
almost  uniformly  along  its  entire  length. 

Radial  profiles  of  the  magnetic  field  components  at 
the  axial  location  indicated  on  Figure  1.  The  axial 
and  radial  components  concentrate  near  the  axis. 

Axial  profiles  of  the  magnetic  field  components  at  the 
radial  location  indicated  in  Figure  1.  The  radial  and 
axial  components  are  largest  near  the  base,  while  the 
azimuthal  component  only  varies  by  a  small  amount 
over  most  of  the  cylinder  length. 

Flux-surface  contours  from  Zweibel  and  Boozer  (1985) . 
The  present  numerical  computations  agree  with  this 
analytic  result. 
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